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Abstract. It is possible to extract, from the observations, distribution functions of the birth 
dynamical properties of a stellar population, and to also infer that these are quite invariant to 
I I ' the physical conditions of star formation. The most famous example is the stellar IMF, and the 
initial binary population (IBP) seems to follow suit. A compact mathematical formulation of 
' the IBP can be derived from the data. It has three broad parts: the IBP of the dominant stellar 
' population (0.08 — 2Mq), the IBP of the more-massive stars and the IBP of brown dwarfs. 
, These three mass regimes correspond to different physical regimes of star formation but not to 
{Dlh' structure in the IMF. With this formulation of the IBP it becomes possible to synthesise the 
stellar-population of whole galaxies. 
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^ [ 1. Introduction 

' The fundamental dynamical properties of stellar populations are the masses of the stars 
and their correlation in multiple stellar systems. The distribution of stellar masses at 
birth, the IMF, is rather well constrained and has (surprisingly) been found to be invari- 
I ant despite theoretical models predicting systematic variation for example of the mean 
stellar mass or even of the minimum mass with the physical conditions of star formation. 
T—i • This problematical issue has been discussed at some length by Kroupa (2008) where the 
IMF Universality Hypothesis is stated. 
, Equivalently, the question may be raised whether the other distribution functions char- 

r^J ' acterising a stellar population, namely the distribution functions of binary systems, are 
■ just as invariant. If this were the case then it would have important bearings on the 
theory of star formation as the fragmentation length-scale may then not depend much 
on the physical conditions of the molecular cloud core. A change in the properties of 
the binary-star distribution functions with mass scale, if found, would yield important 
clues to the fragmentation and angular momentum re-distribution processes during star 
formation. 

The three important distribution functions describing the initial binary population 
(IBP) are the distribution of periods (P, here always in days), or equivalently of semi- 
major axes (a, in AU), the distribution of mass-ratios (q = m2/mi ^ 1) and the dis- 
tribution of orbital eccentricities (e). These are related by Kepler's third law: a^/Pyj. = 
mi -|- m2, where Pyr is the orbital period in years (P = 365.25 Pyi) and rni,m2 are the 
primary- and secondary-star masses in Mq. Because the periods of binary stars range 
over many orders of magnitude the shorthand IP = log^gP is used throughout this text. 
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2. Star formation and the initial binary population (IBP) 

Observations have shown that the star-formation process is intrinsically linked to the 
production of binary stars. Indeed, binary stars must be the dominant formation channel 
because the observed multiplicity fraction (the number of multiple systems divided by 
the number of sources in the survey) is indistinguishably high among old metal-poor 
( Carney et al. (2005) ) and among thin-disk main-sequence stars 

(Duquennoy & Mayor (1991)), and is near unity for pre-main sequence stars and proto- 
stars (Duchene (1999) Connelley et al.(2008) I. If, on the other hand, higher-order multi- 
ple systems were a major outcome of late-type star formation, then the dynamical decay 
of these on a time-scale < 10^ yr would pollute the pre-main-sequence stellar population 
with single stars which are not observed in large numbers. Indeed, this is evidently the 
case for massive stars (m^ few Mq) which appear to form preferentially in binary-rich 
dense cores of populous embedded clusters which rapidly decay dynamically by eject- 
ing massive stars (Clarke & Pringle (1992) Pflamm-Altenburg & Kroupa (2010) I. Thus, 
according to the Binary-Star Conjecture or Theorem ( [Kroupa (2008)[ ) the vast 
number of stars form as binaries, while non-hierarchical higher order multiple systems 
cannot be a significant outcome of late-type star formation. 

The formation of binary systems remains an essentially unsolved problem theoreti- 
cally. Fisher (2004) shows analytically that isolated turbulent cloud cores can produce 
an unquantifiable fraction of binary systems with the very wide range of orbital pe- 
riods as observed. But direct cloud collapse calculations are very limited in predict- 
ing binary-star properties owing to the severe computational difficulties of treating the 
magneto-hydrodynamics together with correct radiation transfer and evolving atomic 
and molecular opacities during collapse. 

The currently most advanced hydrodynamical simulations have been reported by 



Moeckel & Bate (2010) They allow a turbulent SPH cloud to collapse forming a sub- 



stantial cluster of 1253 stars and brown dwarfs amounting to 191 Mq. The cluster has 
a half-mass radius of about 0.05 pc and contains a very substantial binary and higher- 
order multiple stellar population with a large spread of semi-major axes but peaking at 
a few AU. After dynamical evolution with or without expulsion of the residual gas the 
distribution of orbits ends up being quite strongly peaked at a few AU with a signifi- 
cant deficit of orbits with a > 10 AU, and with a deficit of systems with a mass ratio 
q < 0.8, when compared to the main-sequence population (their fig. 11). This state-of-the 
art computation therewith confirms the above stated issue that it remains a significant 
challenge for star-formation theory to account for the Gaussian-type distribution of a 
spanning 10~^ — 10^ AU as for Galactic-field binaries. One essential aspect which is still 
missing from such computations is stellar feedback which starts heating the cloud as 
soon as the first proto stars appear. These heating sources are likely to counter the grav- 
itational collapse such that in reality the extreme densities are not achieved allowing a 
much larger fraction of wide binaries to survive. 

More general theoretical considerations suggest that star-formation in dense clusters 
ought to have a tendency towards a lower binary proportion in warmer molecular clouds 
(i.e. in cluster-forming cores) because of the reduction of available phase-space for binary- 
star formation with increasing temperature (Durisen & Sterzik (1994) hereinafter DS). 
On the other hand, an enhanced binary proportion for orbital periods ZP;$5.6 may be 
expected in dense clusters due to the stimulation of binary formation through tidal shear 
(Horton, Bate & Bonnell (2001)1, thus possibly compensating the DS effect. The initial 
period distribution function (IPF) may thus appear similar in dense and sparse clusters, 
apart from deviations at long periods due to encounters and the cluster tidal field. 
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The multiple-star population in the Galactic field is build-up by star-forming events 
in star-clusters or groups containing from a dozen to possibly millions of stars. Indeed, 
a certain but presently not well known fraction of stars form in small-A'^ systems that 
typically have a size « 100 AU, and the dynamical decay of these is likely to affect the final 
distribution of P^IO^ binaries (Sterzik & Durisen (1998)1, giving rise to non-uniform 
jet activity (Reipurth (2000) I. But again, quantification of their properties is next-to- 
impossible given the neglect of the hydrodynamical component. But, by the Binary- 
Star Conjecture above, the binary formation channel must be vastly dominating 
over the formation of non-hierarchical higher-order multiples. 



3. The IBP universality hypothesis 

It is hoped that star-forming simulations will allow essential insights into reproducing the 
stellar population stemming from an individual modest star- forming event. But given the 
computational complexity, synthesizing for example the binary population in a massive 
star cluster or of a whole galaxy or even parts thereof such as the solar neighbourhood are 
not possible. Therewith it becomes rather apparent that current theory has no predictive 
power concerning the binary properties of stellar populations in different environments. 

However, by proposing that initial distribution functions of the binary-star properties 
(see end of §[T|) exist, that is, that there exists an outcome of the isolated (i.e. low-density, 
such as in Taurus- Auriga) star-formation process that can be quantified in terms of an 
IBP, we would be put into the situation of being able to synthesise populations. This 
would become feasible if it is understood how this IBP is affected by physical processes 
that are inherent to a binary system and that are due to stellar-dynamical encounters in 
denser star- forming regions. 

In fact, the invariance of the IMF must be implying an insensitivity of the star- 
formation outcome to physical conditions. The IMF being invariant constitutes a statisti- 
cal statement on one of the birth dynamical properties of stars (namely their distribution 
of masses). So, since both the IMF and the IBP are the result of the same (star- formation) 
process, and since the IMF is a result of this process "one level deeper down" than the 
IBP, it is quite natural to suggest that the formal mathematical distribution function of 
all of the birth dynamical properties of stars are invariant. Thus the following hypothesis 
follows: 

The Star-Formation Universality Hypothesis: 
IMF universality ■^=> IBP universality. 



4. Inferring the IBP for m < 2 Mq stars 

The challenge of inferring the period-, mass-ratio- and eccentricity-distribution functions 
characterising the IBP can be formulated as follows: assuming the Star- Formation 
Universality Hypothesis to be vahd and using observational constraints on the pre- 
main sequence binary population, the observed main-sequence Galactic-field binary-star 
distribution function must be corrected for the dynamical processes acting in the birth- 
groups or birth-clusters of stars as these emerge into the Galactic field. Fig. [T] visualises 
this idea. 



These dynamical processes are well understood and are detailed in Kroupa (2008) 



Energy arguments imply the Heggie-Hills Law according to which the wide-binaries 
are disrupted preferentially compared to short-period binaries and tight binaries typically 
become tighter. The boundary between the short and long-period binaries depends on 
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the velocity dispersion in the birth population. At the same time, binaries with a small 
mass ratio (small q) are also preferentially disrupted. The tightening of tight binaries 
implies that such systems and individual stars interacting with them can be ejected 
from the cluster. Pre-main-sequence Eigenevolution evolves the short-period binaries 
within a time-scale of few 10^ yr due to system-internal dissipative processes such as tidal- 
circularisation, primary-star-disk-secondary-star interactions and disk-disk interactions. 

The dynamical disruption of binaries induced through the birth cluster and charac- 
terised by the stellar-dynamical operator fj^'^oci.^.s ^ g^j^jj mass-ratio, eccentricity and 
period evolution induced through pre-main sequence eigenevolution, can both be calcu- 



lated (see Kroupa (2008) for details). M^d is the stellar mass of the embedded cluster 
with half- mass radius -R0.5. Note that this "embedded cluster" does not need to be a 
cluster which survives the first few 10 Myr as a bound star cluster. It can readily be 
taken to include sparse star-forming aggregates of stars as observed e.g. in the about 
0.5 pc radius sub-groups of a dozen pre-main sequence stars in Taurus- Auriga. 

The above two-stage transformation of the formal mathematical birth distribution 
functions (/p,/g,/e) to the final distribution function (e.g. of the dispersed cluster) can 
be written as 

2?outcomo(/P, e, (? : mi) = f^A^-i.iio.s [n^.^^^^^^^ [VhiMlP, e, q : mi)]] . (4.1) 

^^birth = fpi^'P) fq(,<l) feie.) is the birth distribution function taking the birth period, 
mass-ratio and eccentricity distribution functions to be separable, that is, the birth pa- 
rameters Pb,qb, Sb are not correlated. 

A simple theoretical treatment of eigenevolution which is based on pre-main sequence 
tidal-circularisation theory has been shown to quite nicely account for the correlations 
between eccentricity, mass-ratio and period for short-period (P ^ 10"^ d) binaries. The 
operator ricigonovoi generates the correlations between period, eccentricity and mass ratio. 

The stellar-dynamical operator, fl^^oci.JSo.s ^ jg given by an Nbody star-cluster, and can 
be envisioned as a transformation of the number of binaries in different binding-energy 
bins such that the reduction is largest in the most weakly-bound binaries. The tightening 
of binaries can increase the fraction of short-period binaries in a population. 

Putting this together, the birth binary population is defined by the following formal 
mathematical rules: 



The birth binary population (BBP): 

• random pairing from the canonical IMF for 0.08 < m/M© < 2; 

• thermal eccentricity distribution of eccentricities, /e(e) = 2e; 

• the period distribution function 

IP ~ IP ■ 

tin mm / . (..N 

fPMrth = V 2 ^ (4-2) 

0+ (IP - iPmin) 

where 77 = 2.5,(5 = 45, ZPmin = 1 and j/jT'""'' fpMrth dlP — 1 such that the birth binary 
fraction is unity (ZPmax = 8.43). 



Details arc provided in Kroupa (2008) Note that customarily the "initial binary popu- 
lation" (IBP) derives from the BBP after pre-main sequence eigenevolution. Given that 
pre-main sequence eigenevolution occurs mostly for short-period binaries on a time-scale 
of 10^ yr within the system, it is expected to be approximately universal. The Star 
Formation Universality Hypothesis thus remains valid. 

Passing the above easily generated birth distributions through eigenevolution and then 
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letting the stellar-dynamical operator J7*^<!oi.^^o.5 g^^^ on this resulting initial distribution 
leads to the Galactic field population as observed. This is shown in Figs [2H31 for the case 
Meci = 13OM0,i?o.5 = 0.8pc, where the arrows indicate the action of r2Afeci,-Ro.5_ These 
results are valid for long-lived cluster models without gas. Short-lived embedded clusters 
would have a smaller i?o.5 to give a dynamically- equivalent J7*^<=<:i'^o.5 (Kroupa (1995a)). 



Note that in Fig. [5] the normalisation of the distribution functions follows the custom 

f=^^, (4.3) 

J Vsys 

where A'orbits is the number of binary-star orbits found in a sample of A^sys — ^orbits + 
.^singles systems or sources, and the A'orbits may be the total number of binaries (for the 
total binary fraction in a population), or just the number of orbits in an IP interval 
yielding fp. 



5. The IBP for massive stars (m > 2Mq) and for brown dwarfs 

A seminal argument suggesting that massive stars form in regions of high density void 
of low-mass stars but preferentially in tight binaries with similar component masses has 



been provided by Clarke & Pringle (1992) They reach this conclusion on the basis of 



the distribution of OB runaway stars finding these conditions to be necessary. In Bonn, 
Seungkyung Oh is building on these semi-analytical results by performing direct Nbody 
calculations of realistic young clusters and testing different pairing rules for massive stars 
in initially mass-segregated and unsegregated clusters. The general result is that massive 
stars indeed need high densities and tight binaries with mass-ratios nearby unity in order 
to account for their spatial distribution around star forming regions. 

Brown dwarfs (BDs) form an altogether distinct population from stars as is deduced by 
Thies & Kroupa (2008) on the basis of the observed different pairing properties of brown 
dwarfs and very-low-mass stars on the one hand side, and stars on the other hand side. 
The most famous evidence for this comes from the brown-dwarf desert: while stars pair 
up irrespectively of their masses (e.g. G-dwarfs typically have M-dwarf companions), BDs 
are exempt from participating (there are exceptions, but such exceptions are a natural 
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Figure 1. Schematic of how the observed Galactic-field main-sequence (thick black curves, 
Duquennoy & Mayor (1991) I period or binding energy distribution (left panel) and the 
main-sequence mass-ratio distribution (right panel) are back-computed (upwards pointing ar- 
rows) with the stellar-dynamical operator (1*^"'=''^" to yield an estimate of the initial (thin 
red curves) period and mass-ratio distributions. Note that the physical reality of Q^'^^t^i'^o.s jg 
established if both fp and fq become consistent with the pre-main sequence data for one and 
the same f^^^od.flo.s 'pjjg eccentricity distribution is not afi'ected by r2*^<=<:i.«o.5 
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outcome of stellar-dynamical processes). BDs also have a very different semi-major axis 
distribution, by being limited to a ^ 20 AU, while for stars a ranges over five orders of 
magnitude. Also, BDs have a small binary fraction / « 0.15 with a bias towards similar 
companion masses. Computing the effects of the pairing properties of BDs on the observed 
mass function, Thies & Kroupa (2008) deduce that the IMF must be discontinuous near 
0.08 M©. 

Massive stars and BDs thus cannot be included in the BBP formalism above, but 
require their own mathematical rules. In terms of physics this means that the formation 
of massive stars and brown dwarfs occur in a physical regime different to that of the 
typical star. 

This is naturally understandable by noting that massive stars can only form in the 
densest regions of their parent embedded cluster. BDs, on the other hand, typically form 
in the outer regions of extended circum-stellar discs of the average star, either by disc 
instability (Goodwin & Whitworth (2007) ) or by tidally-induced gravitational instability 
caused by passing stars in the parent cluster (Thies et al. (2010) I, or are ejected embryos 
( [Reipurth fc Clarke (20(5I)| ). 

In terms of numbers, massive stars and brown dwarfs are much rarer than the typical 
star: on average one BD forms per five late-type stars, while one massive star forms per 
few hundred late-type stars. 
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Figure 2. Left panel: The transformation of the birth period distributio n function (BPF, 
Eq. [42I thick red curve with = 2.5,6 = 45, ZPmin = 1, 

eigenevolved "initial" period distribution function (IPF, dasE 

described by Eq. lU with « 3.5, 5 « 100, IPmin = 0, [Kroupa (1995a) 
histogram is the final PF after ^isoMg.o.Spc ^^^^ jgp^ rpj^^ ^^jj^ ^^^^^ ^p^^ circles and 

stars are G-, K-, an d M-dwarf binaries, and the open squares are pre-main sequence systems 
(see [Kroupa (2008) | for references). The normalisation of this plot is such that each period 
bin contains the fraction of binary orbits in the whole sample of stars plus binaries. Right 
panel: The transformation of the initial mass-ratio distribution for primary masses mi « 1 Mq 
(dashed histogram) to the final mass-ratio distribution after the population emerges from its 
star cluster (solid histogram, from K2). The upper panel is for short-period binaries which are 
not affected by fi''-^" Mq,o.Spc ^ while the bottom panel shows the long-period binaries. Note that 
the birth mass-ratio distribution, which results from random pairing from the IMF, is given 
by the dashed histogram in the lower panel, while eigenevolution transforms this distribution 
to the dashed histogram evident in the upper panel. That pre-main sequence, i.e. dynamically 
unevolved, mass-rati os are consistent with random paring from the IMF has been found by 
[Woitas et al. (2001) [ 
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6. IBP mass regimes and IMF structure 

A possibly interesting issue thus emerges: While the BD/star schism is clearly evident 
in the different pairing rules and the discontinuity in the IMF near 0.08 M0, the change 
of pairing rules between late- type and early- type stars (roughly at a few Mq) does not 
seem to be evident in any corresponding structure in the IMF. In fact, the flattening 
of the IMF near 0.5 M© does not seem to correspond to a change in birth binary-star 
properties. While this fact does not contradict the above statements, it does provide an 
additional constraint on star-formation theories. 



7. Dynamical Population Synthesis (DyPoS) 

Having thus obtained a formal mathematical description of the invariant birth binary 
population it now becomes possible to perform Dynamical Population Synthesis 
(DyPoS). DyPoS rests on the same ansatz as has already been applied to model the 
vertical structure of the Milky Way disk (Kroupa (2002) ) and to model the stellar initial 
mass function of a whole galaxy (Kroupa & Weidner (2003) ) by adding up the contribu- 
tions by each embedded cluster. This Lego Principle comes about from the realisation 
that stars that ultimately end up populating the field form in groups embedded in gas 
(e.g. |Lada fc Lada (2003)[ ). 

If the birth stellar population in the groups and the transformations that occur within 
the groups before the population hatches into the field are understood and are known, 
then an entire galaxy can be synthesised. This principle is easy to implement for the inte- 
grated initial mass function (IGIMF) . It is a little bit more involved when implementing 
it for an integrated description of the vertical structure of the Milky Way disk because 
the stellar velocity field emerging from an embedded stellar group needs additional as- 
sumptions on its time-evolving properties. 

Along the same lines, the binary population in a galaxy follows from adding up all 
the binary populations born in the embedded groups taking into account the dynamical 




q eccetitridity 



Figure 3. Left paneh The transformation of the overall initial mass-ratio distribution for 
0.8 < mi/MQ < 1 (dashed histogram) to the final mass- ratio distributio n after the popu- 
lation emerges from its star cluster (solid histogram, from [Kroupa (2008) | ). The q = 1 peak 
results from pre-ma in sequence eigenevolution. The solid circles are observational data by 
|Reid fc Gizis (1997)[ Right paneh The eccentricity distribution for short-period {IP < 3, 
dashed histogram) and long period {IP > 3, solid histogram) binaries. The thick red solid curve 
visualises the bell-shaped distribution of short-period orbits which results from the thermal dis- 
tribution (the solid histogram) after eigenevolution. Stellar-dynamical encounters, i.e. fi^^cd'^o-s 
have no effect on the thermal eccentricity distribution (it is an invariant to Sl^^cd-^o s ■) 
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transformation that acts on the population in each group as it evolves and disperses, 
Vfi,id{lP, e,q:mi)= / j^A^d.flo.s [J7,ig,„,,„i [I?birth(ZP, e, q : toi)]] 

X ^ecl(-Wecl,-Ro.5) ^^0.5 dM^^h 

(7.1) 

where ^eci(Afeci, -R0.5) is the distribution of stellar masses and half- mass radii of embed- 
ded aggregates of stars ("embedded clusters") forming in a time interval St (« lOMyr) 
throughout the galaxy. Eq. 17.11 sums up all the evolved populations of binaries that are 
formed in ^cci d.Ro.5 dMcd clusters. i?o.5 ~ 0.4 pc for embedded clusters and ^cci oc M~J 
is the power-law mass function (/3 = 2 is the usually found index from observational 
surveys). 

Two examples of DyPoS have already been computed by Michael Marks at Bonn 
University. Eq. 17.11 is solved on a grid of i?o.5 and power-law index, /3, of the embed- 
ded group or cluster initial mass function, ^cci(-^^cci), where the maximal star-cluster 
mass, Meci,max follows from the star-formation-rate (SFR) versus Mc, 



Weidner et al. (2004) In an elliptical galaxy which formed in a burst with SFR 



lO*M0/yr the most-massive star-cluster weighs about 10 Mq corresponding to the 
mass-scale of ultra-compact dwarf galaxies. Binary systems are disrupted efficiently in 
the massive clusters. Using Eg. 17.11 and assuming the typical embedded star cluster has 
a radius of 0.4 pc, the binary fraction of the late-type stellar population is computed 
to be /bin — 0.45. In a dirr galaxy with a SFR = lO~*M0/yr the most massive clus- 
ter that can form has a mass of about AfocUmax — 100 Mq. For R0.5 = 0.4 pc DyPoS 
yields /bin — 0.85, because many more wider binaries survive the on average less-massive 
embedded clusters in the dIrr galaxy, compared to the above star-burst E galaxy. 



A full mathematical treatment is found in Marks, Oh & Kroupa (2010) 



8. Conclusions 

It appears that the star-formation outcome in terms of stellar masses and multiple sys- 
tems can be formulated by the Star Formation Universality Hypothesis (§ 111)- 
For stars with m^2MQ the Birth Binary Population (§ g]) can be defined. This 
is the outcome of star formation in low to intermediate density {p^ IO^'Mq/pc^) cloud 
regions (e.g. of an embedded cluster). For m<^2MQ stars the pairing rules change (§[5]) 
perhaps reflecting the outcome of star formation in dense regions such as in the cores 
of embedded clusters (p^ 10^ MqPC^). Brown dwarfs follow entirely separate rules (§[5]) 
being an accompanying but distinct population to stars. It remains to be understood 
why these changing IBP properties do not correspond to the structure evident in the 
IMF (§111). 
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